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ABSTRACT 


The  scattering  of  radar  energy  from  a rough  ocean  surface 
is  considered  using  the  approximations  of  physical  optics  and 
perturbation  theory.  The  ocean  surface  is  taken  to  be  generated 
by  a random  process  with  zero  mean  value  and  a homogeneous 
covariance  function  represented  by  a wavenumber  energy  spectrum. 

For  slightly  rough  seas  the  coherent  and  incoherent  backscatter 
cross-sections  are  obtained  for  a general  ocean  spectrum.  For 
very  rough  seas,  the  physical  optics  cross-sections  were  calcu- 

i 

lated  for  a Gaussian  ocean  spectrum  neglecting  shadowing.  ! 

! 

Scattering  from  ripples  where  the  ocean  surface  radius  of  cur- 
vature is  no  longer  large  compared  to  radar  wavelength  was 
treated  by  perturbation  theory  using  the  magnetic  dyadic 
Green's  function  for  a half  plane.  The  scattered  field  per 
turbed  by  the  rough  surface  was  obtained  from  a perturbation 
expansion  of  the  surface  fluctuation.  Incoherent  cross-sections  ♦ 

i 

were  found  to  be  different  for  horizontal  and  vertical  polar- 
izations. The  "visibility"  of  current  induced  perturbations 
for  particular  radar  frequencies  is  calculated  for  a Phillip's 
saturated  spectrum  of  surface  waves. 
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I . INTRODUCTION 


The  standard  treatments  of  radar  scattering  from 
the  ocean  surface  generally  assume  i)  that  this  surface 
wave  field  has  an  incoherent  hydrodynamic  interaction  and 
is  therefore  describable  by  a spectrum,  and  ii)  that  this 
spectrum  may  be  partitioned  into  a part  with  waves  of  low 
curvature  and  a part  with  waves  of  low  amplitudes.  Callen 
and  Dashen^  argue  this  separation  at  a wavenumber  kc  by 

comparing  the  mean  square  radius  oi  /ature  for  the  long 

2 2 
waves  R and  the  mean  square  height  o*.  the  diort  waves  h 

to  the  radar  wavelength  Xr.  The  conditions  h/Xr  <<  1 and 

X^R  <<  1 are  then  satisfied  simultaneously.  The  short 

waves  are  the  scattering  centers  and  ride  on  the  longer 

waves.  The  scattering  of  scalar  electromagnetic  waves  from 

the  short  water  waves  is  treated  in  linear  perturbation 

theory  and  the  undulating  surface  on  which  the  scatterers 

ride  by  geometric  optics. 

2 

Calculations  by  Lane  indicate  that  the  dominant 
effect  of  the  longer  waves  is  to  provide  a tilted  surface 

. 3 

on  which  the  scattering  centers  ride.  Wright  calculates 
the  effect  of  the  long  waves  (k  < kc)  by  averaging  over  the 
tilt  angle  of  the  surface  in  determining  the  scattered 
radiation.  The  general  expression  for  the  radar  cross 
section  of  the  surface  can  be  written  as 


I 


the  angle  of  incidence  of  the  radar  (9). 

A second  effect  of  the  long  waves  and  one  included 
in  the  second  term  under  the  integral  is  the  phenomenon  of 
shadowing.  Shadowing  is  a function  of  surface  slope 
rather  than  curvature  and  therefore  has  contributions 
from  both  the  long  and  short  waves.  This  phenomenon  is 
the  blocking  of  electromagnetic  waves  by  one  segment  of 
ocean  surface  so  that  they  do  not  reach  a second  seg- 
ment of  ocean.  If  this  effect  is  dependent  only  on  the 
angle  of  incidence  and  the  local  slope  of  the  surface,  then 
it  may  be  described  by  a function  S(0r$,«f>).  The  agruments  used 
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the  ocean  will  be  in  the  shadow  cf  another  for  a given 
angle  of  incidence. 

In  Beckmann^  this  probability  is  related  to  the 
shadow  function  S(6)  for  scattering  from  a random  rough 
surface.  This  function  is  zero  on  the  shadowed  parts  of 
the  surface  and  unity  on  the  illuminated  sections.  A 
section  of  the  surface  with  unit  area  will  present  a 
reduced  area  S(8)  (<1)  available  for  scattering  when  shadow- 
ing is  present. 

Bass,  et  al.®,  discuss  the  radar  cross  section  as 
a function  of  angle  of  incidence  implementing  various 
assumptions  about  the  statistics  of  the  longer  waves,  i.e., 
k < kc*  A shadowing  function  is  obtained  by  averaging 
over  the  fluctuations  of  the  longer  waves  and  the  sensi- 
tivity to  the  statistical  character  of  the  longer  waves 
is  studied.  The  radar  cross  section  (a)  is  expressed  in 
terms  of  a scattering  integral  which  includes  the  shadowing 
function  in  the  integrand.  This  result  is  formally  exact, 
subject  to  the  general  restrictions  noted  in  the  Callen- 
Dahsen  analysis. 

The  third  term  under  the  integral  in  Eq.  (1)  is  a 
geometric  factor  determined  by  the  boundary  conditions  for 
electromagnetic  scattering  from  an  inclined  surface.  This 
factor  is  clearly  a function  of  the  incident  ray  angle  and 
the  tilt  of  the  surface,  i.e.,  the  slope. 


cm  P^rir 


The  last  factor  in  Eq.  (1)  is  the  spectral  density 
function  for  the  small  wave  structure  at  the  ocean  surface 
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(k  > k ).  It  is,  of  course,  this  structure  which  gives 

%*« 

rise  to  the  backscattered  radiation  and  provides  the 
"random  rough  surface"  assumed  to  be  present  in  the  dis- 
cussions of  the  shadowing  function.  These  explicit 
assumptions  abouc  the  independent  and  stationary  character 
of  the  statistics  of  the  short  and  long  waves  are  question- 
able. The  modulation  of  the  short  wavelength  spectrum  by 

waves  near  the  spectral  peak  is  discussed  in  Watson  and 
7 

West  , but  no  systematic  analysis  of  the  effect  on  the 

statistics  of  the  short  waves  yet  exists. 

The  treatments  of  shadowing  mentioned  above  have 

2 3 5 6 

all  neglected  diffraction  ' ' ' ' The  effect  of  diffrac- 
cion,  using  the  mean  square  slope  of  the  ocean  surface  as 
a parameter,  is  discussed  in  Appendix  A. 


In  this  report  we  shall  be  concerned  in  part  with  the 
vector  scattering  of  radar  waves  from  a random  ocean  using  the 
approximations  of  physical  optics  and  perturbation  theory. 

The  surface  C(x,t)  is  taken  to  be  a random  function  continuous 
in  the  mean  with  a specified  correlation  function  of  surface 
heights.  The  correlation  function  is  characterized  by  the  sea 
surface  spectrum  for  the  short  wavelength  waves.  The  scattered 
far  field  is  related  to  the  incident  radar  wave  by  an  approx- 
imate form  of  the  Stratton-Chu  theory. 
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In  the  physical  optics  approximation,  the  basic  assumption 
is  that  every  point  of  the  rough  ocean  surface  reflects 
the  incident  wave  like  an  infinite  tangent  plane.  At 
radar  frequencies,  the  surface  is  perfectly  conducting. 
Shadowing  and  multiple  scattering  effects  are  neglected. 


5 


II.  PHYSICAL  OPTICS  SCATTERING 

The  electromagnetic  field  scattered  from  an  area 
of  ocean  S is  given  by8 

?(r)  = J dS  |nx|v'  x[^Q(r' ) + f (?•  )]j  •'fel?-?' ) 


(2.1) 


where  i£(r)  is  the  scattered  field  at  receiver  point  r,  S 

A. 

the  illuminated  ocean  surface,  n the  local  surface  normal,  and 
£Q(r')  the  incident  radar  field  at  a point  r'  on  the  surface. 
The  free-space  dyadic  Green's  function  is 


G(r-r' ) 


= ( ^ + Tf)  exp(-k|r-r*  |)/4tt  |r-r'| 
v k * 

r-»-«>  (*<f  - rr  J exp£-ik  (r-r  «r' )]  /4irr 


(2.2) 


The  scattering  geometry  is  shown  in  Figure  1. 

If  the  far  field  approximation  is  used  in  Eq.  (2.1) 


we  obtain 


Ei(?)  = 2Frexp(_ikr)  kEoPij  / dS  "j  exp  [-i  (Jt-kr)  .r '] 


pij  - «« 


- : ; ) c £(o)  . ; s(o) 

ij  rirj’  lkiEj  kjEi  J 


(2.3) 


where  the  tensor  P^  contains  all  polarization  information. 
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We  note  that  due  to  the  partitioning  of  effects,  ensemble 
averaging  can  be  applied  directly  to  the  surface  integral 
in  Eq. (2.3)  free  of  polarization. 

To  express  the  surface  integral  £call  it  ?(r)j 
in  terms  of  the  surface  height  random  variable  r(x'),  the 
surface  displacement  is  assumed  to  be  a first  order  quantity, 
i.e.,  the  long  waves  are  essentially  flat  so  that 


r'  = xx ' + yy ' + z^(x') 


n = V'lz'-Cfx'  )j 
| V(z'-£) | 


(2.4) 


and  an  element  of  area  on  the  inclined  surface  is 


dS  = ( V(2'-5) | dx'dy' 


(2.5) 


When  Eqs.  (2.4)  and  (2. 5) are  used  in  the  stochastic  integral, 

J(r)  = fdx’dy’  ( z-x  - y ~^rj  exp£-i  (ic-kr)  . [xx'+Jy’+z^  (x*  )]| 


(2.6) 


For  the  important  case  of  backscatter,  jc-kr  = 2Jc,  and  de- 

/V  A 

noting  the  incident  polarization  by  EQ  = q,  we  obtain 
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III.  PHYSICAL  OPTICS  SCATTERING  FOR  SLIGHTLY  ROUGH  SEA 

The  integral  in  Eq.  (2.7),  call  it  I,  is  the  basis 
for  our  study  of  radar  scattering  from  a rough  ocean.  If 
we  square  and  average  the  scattered  field,  we  obtain  the 
mean  backscattered  power  as  a function  of  incidence  angle 
0Q.  We  note  that  correlations  of  heights  and  slopes  of 
the  sea  surface  must  be  known  to  evaluate  the  backscattered 
power.  In  this  section  we  assume  the  roughness  height  to 
be  small  (ko  <<  1)  and  any  correlation  distance  to  be 
much  greater  than  the  wavelengths  of  the  scattering  centers. 

The  mean  radar  cress  section  is  defined  by 


(3.1) 


and  can  be  expressed  as  the  sum  of  the  coherent  (a  ) and 

c 

incoherent  (a-j.)  cross  sections 


(3.2) 


Averaging  over  the  square  of  the  integral  part  of 
Eq.  (2.7)  yields 
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<lll2>  = / exp^i2k  sinQg  (x^-x2)j  dx-^cb^dy^y. 


cos20o  +[^21)-^22)]  cos6o  sin6o  +^3  sin20o 


(3.3) 


If  the  statistical  fluctuations  in  the  surface  displacement 

^5(x)J  are  homogeneous  and  Gaussian  with  zero  mean,  then 

2 

in  terms  of  the  correlation  function  oz  Efx-^-x^,  y-^-y2)  / 
we  have 

£ = ^exp  J^i2kcos0o(C1“C2)]^=  exp p(2kazcos©o)  2 (1-E ) j 


^21]  E exp[i2kGOS0o(crc2)]^ 

i 

x exp  (2kazcos8o)  2 (l-E)j 
63  5 ^ exP[  i2Kcos0o(?1-C2!]>  = 


i2ka 


2 3_E_ 
z 3x. 


,2  3' 


z 3x^3x2 


+ (2kcos0 


o°z)2  "z  exp[-(2kci2coS0o)2a-E)] 


(3.4) 


where  is  obtained  by  interchanging  1 and  2 in  the  expres- 
sion for  £>2^  and  Sj,  = 1 = 1'2* 

Using  these  results  in  Eq. (3.3)  we  obtain  after 
some  reduction 
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<IH2>  = exp|-(2kazcoseo)2]  J <3xidx2dy1dy2  exp£i2ksineo(Xj-x2)J 


.2.  . !!li°  _Z_  * jl 


- eQ  - + — 


exp  ({2kaz 


coseo) ‘E 


(3.5) 


A similar  computation  shows 


<I>=  cos0Q  exp[-)s(2  o2coseo)2]  / dx-^dy^  exp(i2  sineox-^)  (3.6) 


where  we  have  used 


exp(i2kccos0Q)^  = exp  -*s  (2kazcos0Q)  2] 
exp(i2kCcos0Q)^  = 0 


For  small  root  mean  square  (rms)  fluctuations  in  surface  height. 


a = (2kozcos0o)  <<  1 , 


Eq.  (3. 5) can  be  approximated  by 


^|I|2^>  a exp  (-a  )j  dx1...dy2  ^ 


2 

cos20  (1+aB  - (o  sine  )2 

0 z ° 3x^2 


2 2 
+ i4ka  sinO^cos  0 
z o 


« 

o ^r]  Okp[i2ksineo(x1-x2)j 


.12 


(3.7) 
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The  two-point  correlation  function  E can  be  expressed 
as  the  Fourier  Transform  of  the  energy  density  of  ocean 
surface  waves  <5>(k)  , 

CD 

E(x1-x2,y1-y2)  = J d <$(<)  exp^i  (x^x^ +<2  (y^y^jj  . (3.8) 

— 00 

Using  this  representation  and  assuming  the  illuminated 
surface  area  to  be  large  in  both  dimensions,  we  obtain 
after  integration  over  S 


^|l|2^exp(a)  = cos20  C2(x2k  sine  ) 


+ acos2eo  I d<  <j>(ic)  C2(ic+x2k  sineo) 


+ (ozsin0 


CO 

J" d k <Mic) 

— CO 

00 

%> 2 J a*  4 


4>(<)  C ( x+x2k  sin0o) 


’°/d 


- 4k  a2  sin8Q  cos20q  / d<  <J>(k)  C2(<+x2k  sin0  ) 


where 


:(ic)  = J dxdy  exp  Ji  (k-jX+k^)] 


(3.9) 


We  note  that  C(k)  is  sharply  peaked  about  tc  = 0 for  a large 
surface  and  approaches  a delta  function  as  S + ®.  Making 
use  of  this  property,  we  obtain  the  following  approximation 
for  Eg.  (3.S)  : 


^|I  |2^>exp(a)  = cos2eo  C2(x2k  sin0Q)  + (2tt)  2 S (2kaz)2  $(~x2k  sin0Q) 


(3.10) 


On  squaring  Eq.  (3.6)  we  obtain 


KOI2  = cos20o  exp  (-a)  C2(x2k  sin0Q) 


(3.11) 


so  that  using  Eqs.  (3. 10)  and  (3.11)in  Eq.  (3.2),  the  incoherent 
cross  section  per  unit  area  is 


-<?■  = Ion  (k2oz)2  exp£-(2kazcos0Q) 2 J $(-x2k  sin0Q)  . (3.12) 


The  coherent  cross  section  per  unit  area  is,  similarly. 


-f  = IT  cos2qo  exp[-(2kazcos0o)2J 


C (x2ksin0o) 


(3.13) 


The  rough  ocean  surface  is  much  larger  than  the  area 
illuminated  by  the  radar  antenna  beam.  The  field  strength 
over  the  illuminated  area  is  not  uniform  but  varies  according 
to  the  gain  pattern.  We  shall  assume  that  the  field  is 
uniform  over  the  half  power  angular  beam  width  0_  and  zero 
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outside.  For  moderately  large  angles  of  incidence,  the 
geometry  is  shown  in  Figure  2.  The  radar  pulse  is  t 
seconds  in  duration,  c is  the  velocity  of  light,  and  r is 
slant  range  to  the  illuminated  area.  The  rough  surface 
area  contributing  to  the  backscattered  power  is  one-half 
of  the  illuminated  area  or 


S = h Aj  = % c6Brx 


If  the  illuminated  area  is  approximately  rectangular  with 

CT 

dimensions  xQ  = r0B  and  yQ  = -y  , then 


C(x2k  sin0Q)  = S sinc(kxQ  sin0Q) 


and  the  coherent  cross  section  per  unit  area  is 


~ - cos20Q  exp  (2kazcos0o)  2 J sine2 (kxQ  sin0Q) 


Since  Xr/xQ  <<:  1»  the  coherent  cross  section  is  unimportant 
off  normal  incidence. 


- 


CT 


FIG.  2 

RmDAR  pulse  illumination  of  ocean  surface 
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IV.  COMPOSITE  MODEL  - SCALAR  EQUATION 


In  Section  II,  it  was  shown  that  the  effect  of  polari- 
zation factors  from  the  integral  deteirmining  the  scattered 
electric  field  [see  Eq.  (2.7)].  A direct  evaluation  of  the 
scattering  integral,  which  is  a scalar,  in  terms  of  an  en- 
semble average  was  than  made  in  Section  III  to  determine  the 
dependence  of  the  scattered  field  on  the  spectrum  of  scatterers 
and  the  undulating  surface.  In  this  section,  we  extend  the 
investigation  of  Callen  and  Dashen^"  on  the  scattering  of  a 
scalar  wave  from  a rough  surface  to  obtain  an  explicit  de- 
pendence of  the  clutter  cross  section  on  the  spectrum  of 
long  waves.  The  analysis  follows  that  cf  reference  (1)  until 
the  evaluation  of  the  scattering  integral  is  made.  At  this 
point,  Callen  and  Dashen  use  the  stationary  phase  approxi- 
mation to  evaluate  the  integral,  whereas  we  make  explicit 
assumptions  about  the  statistics  of  the  long  wave  amplitudes 
and  about  the  energy  spectrum  to  simplify  the  integral. 

In  the  investigation  of  Callen  and  Dashen  the  vector 
field  given  by  Eq.  (2.1)  is  replaced  by  the  scalar  equation, 

E(1)  (x)  = -k2(n-l)y*dp  CX(P)  Go[L?(P)]E{o)[x(p)]  (4.1) 
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where  the  total  field  is  given  by  the  perturbation  expansion 


<?>  =E 


E(P>  (?) 


(4.2) 


E ? being  the  surface  displacement  and  n the  dielectric 

constant  for  sea  water.  The  surface  has  been  separated  into 
a long  and  short  wavelength  part  with  displacements  ?q(p)  and 
C-^(p),  respectively,  as  shewn  in  Figure  3.  The  vector  p = (x,y) 
locates  a point  on  the  surface  z = 0,  and  x(p)  - (P/Cq(p)) 
locates  a point  on  the  undulating  surface  z = CQ(p}. 

Using  the  Greens'  function  from  reference  (1)  and  the 
unperturbed  scalar  field  in  Eq.(4.1),  we  obtain  the  far  field 
expression. 


EU)(k'R)  = J dp  Ci 


with  the  geometric  factor 


(p)  T (a)  T(a')  exp  Jik (k-k' ) * x(P)J 


(4.3) 


T(a)  = 


2 cosa 


cosa  + / .2 

^n-sin  a 


(4.4) 


arising  from  the  boundary  condition  for  the  electromagnetic 

A A 

wave  at  the  undulating  surface.  For  backscatter  k'  = -k 

A A 

and  cosa'  = cosa  = k*  *n,  so  that 


T(a)T(a')  ~ 


4 cos  a 
n + 2/n~ 


where  the  denominator  differs  from  the  exact  expression  by 
less  than  1%  for  n = 80. 

The  integral  inEq.  (4.3)  may  be  changed  to  a surface 
integral  over  the  illuminated  surface  CQ(p)  by  introducing 
which  is  the  normal  distance  between  the  surface 
and  + z^.  We  then  have  dp?1(p)  -*•  dS  ^(x^y).  Letting 
Z ( x , y ) = ZQ(x,y) , we  have  for  the  unit  normal  to  Z(x,y), 


z - x Z - y Z 
x 1 \ 


2 2 

1 + Z + Z 
x y 


the  unit  area 


dS  = / 1 + Z^  + Z^  dxdy 
x y J 


the  cosine  of  the  local  angle  of  incidence 


~ Z„  sin0_  - cos0 


cosa  = k • n = 


2 2 
1 + Z + Z* 
x y 


;X(P)  = ?1(P)  / 1 + Z2  + Z2  . (4.9) 


Using  Eq. (4.7)  - (4.9)  in  eq  (4.3)  yields  for  the  back- 
scattered  perturbed  field 


Et;L>  ,-KB)  - *t<-n-» 


n (n+2/n) R 


J dy  ^(X/y)  £cos20q  + Z2  sin20c 


- Zx  sin20Qjx  exp  £-i2k(x  sin0Q+Z(x,y)  cos0Q)j . (4.10) 


£1 


I I 


Using  Eq.  (4.10),  we  introduce  the  function  fQ(p)  defined 


by 


f (p)  5 (cos20  + Z2  sin20  - Z sin20  ) exp -i2k(x  sin0  + Z cos0  ) 

O OX  OX  O L O OJ 

(4.11) 

so  that  upon  squaring  and  averaging  (4.10)  over  an  appropriate 
ensemble  we  obtain 

<|E‘1,|2>  = AteT/a?/  d5,<Cl,5,tl,J,,><*o‘»>£o,5'»>  • 


(4.12) 


In  Eq.(4.12)  we  have  assumed  that  the  ripples  (5  ) are 
statistically  independent  of  the  long  waves  ( f ) . 

In  the  following  analysis,  we  assume  a Gaussian  ocean 
surface  for  the  long  waves  and  that 
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V<p>  Z(p')>=  cr  Ez(p-^) 

where 

°*2  * <z2<*>> 

and  Ez  is  the  correlation  function  for  the  long  waves. 
Using  the  notation  z = z{p>,  z'  - z(p')  and  K = -2k  cos  0 

c 

we  have 

<fo«»fS  exp[-i2k(x-x')  sinej  j cos40o  ^expjiK(z-Z') 

+ sin2eocos20o  <{z2  + z'?)  exp[iK(z-z')J) 
cos  0Q  Sin20o^(Zx  + Z^,)  exp|iK(Z-z')J  ’ 
sm  0q  sm2©o  ^ZXZ^/ Ux'+  zx)  exp  jiK  (Z-: 
+ sm40Q  <Z2  2 x,  exp JiK  (Z— Z ' )J  ^ 


+ sin220  z*' 


f ■>  ir(’7  —’7  aO 


The  averages  defined  by  eq  (4.15)  are  evaluated  in  Appendix  J 

i 

B for  gentle  long  waves.  j 

We  introduce  the  function  F(p-p  >m,QQ)  by  the  right  hand  | 

S 

side  of  eq  (4.15)  , as  follows  f 

3 

n 
5 

i 

^o(p)f*{p,)^>  H F(p-p';G>0)  exp  £-i2k  (x-x*)  sin0QJ  (4.16)  j 

| 

Equation  (4.16)  may  be  used  to  evaluate  the  integral  in  (4.12)  ‘ 

3 
s 

when  the  ripple  correlation  length  is  small  compared  to  the  1 

I 

long  wave  correlation  length , i . e . , the  correlation  of  short 
waves  is  essentially  a delta  function  on  the  scale  of  the 
long  wave,  so  that  ; 

<lE<1)l2>'*#r-  F(I;0O)  /a? 

7T  R 

x exp£-i2k  sin0Q(x-x/)j  (4.17) 

where  | l | is  the  correlation  length  of  the  short  waves . The 
integral  defines  the  spectrum  of  short  waves  at  X&  = (2k  sin0Q,  0)  , 
r=p-p' , yielding  for  the  clutter  cross  section 

a(D  = r2<Je(1)  )2> 

s 

I 

I * 2.61k4  FU;0o)  <5>(-Ak,r)  . (4.18) 

[ 

where  S is  the  illuminated  area.  The  cross  section  for 
scattering  from  the  long  waves  as  given  in  reference ^ is 
unchanged. 
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V.  OCEAN  RADAR  BACKSCATTERING  CALCULATIONS 

In  the  preceding  sections  and  in  Appendices  C,  D,  and 
E,  a number  of  distinct  models  for  the  scattering  of  electro- 
magnetic waves  from  ocean  surfaces  of  varying  degrees  of 
roughness  are  presented.  The  models  have  been  based  on  dif- 
ferent combinations  of  the  four  conditions:  (i)  that  the 

roughness  height  be  smaller  than  the  radar  wavelength  illumin- 
ating the  surface  (kaz  <<  1) , (ii)  any  correlation  distance  for 

the  short  waves  is  smaller  than  their  wavelengths  (kJt  <<  1) , 

2 

(iii)  the  local  mean  square  slope  be  small  « 1)  and 

SGcl 

(iv)  the  root  mean  square  radius  of  curvature  of  the  long 
waves  is  much  greater  than  the  radar  wavelength  (Ar/R  <<  1) . 
Condition  (iv)  is  a restriction  in  all  the  models  discussed. 

The  applicability  of  the  above  models  to  various  wave- 
number  regimes  and  sea  states  will  be  explored  in  this  section 
Also,  the  statistical  parameters  which  characterize  the  ocean 
surface  will  be  discussed  for  a variety  of  representations  of 
the  ocean  wave  spectrum. 

In  Figure  4 , we  show  a schematic  of  an  ocean  energy 
spectrum.  We  note  especially  that  the  frequency  region  of 
interest  for  radar  returns  spans  both  gravity  and  capillary 

. 9 

waves.  The  Pierson-Moskowitz  spectrum  gives  a good  represen- 
tation of  observed  properties  of  wind  generated  gravity  waves 
in  deep  water. 


where  in  cgs  units,  g = 980  cm/sec  , U is  the  wind  speed 

(cm/sec)  and  a and  8 are  dimensionless  constants  with  the  values 
-3 

8.1x10  and  .74,  respectively.  Descriptions  such  as  Eq.  (5.1) 
are  not  as  useful  as  they  could  be,  however,  since  most 
theoretical  discussion  of  the  evolution  of  the  surface  wave 
field  is  in  terms  of  the  wavenumber  spectrum'*'0''1'^''1'2.  Using 
the  relation 

<Moj)  du  = medic  («)  (5.2) 

(2) 

where  $>  (k)  is  the  two-dimensional  wavenumber  spectrum 

and  the  linear  dispersion  relation  for  gravity  waves 
o»  = v^jk,  we  obtain 

$f2)  (fc)  = _o_  exp(-eg2/U4<2)  (5.3) 

2ttk 

which,  neglecting  the  exponential,  yields  the  Phillips 

13 

saturated  spectrum 
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In  Project  SWOP  a stereo  technique  from  two  air- 
planes was  used  to  measure  the  two-dimensional  wavenumber 
spectrum.  The  measurements  were  made  over  the  Atlantic 
ocean  (25  Oct  54,  39°N,  63.5°W)  over  an  area  2700  ftxl800  ft. 
Wind  speeds  ranged  over  13-*20  knots  during  the  preceding 
14  hours  so  that  the  spectral  cut-off  frequency  in  Eq.  (5.1) 


ranged  over  .444  to  .10  cps.  An  average  frequency  of  .29 

cps  was  used  to  determine  a cut-off  wavelength  X = = 

or 

Figure  5 shows  the  final  smoothed  spectrum.  Figure  6 shows 

the  one-dimensional  wavenumber  spectrum  taken  along  the 

north-south  direction  in  Figure  5.  Numerical  integration 

of  the  latter  curve  leads  to  a mean  square  surface  height 

fluctuation  of  o =6.75  ft  and  a lower  bound  for  the  mean 
z 

sauare  slope  of  (|)2  _ = 8.35x10  3 . For  < > .02  ft  \ 

" S6a 

-4 

the  data  indicated  an  excellent  < fit.  The  analytical 
expression  for  the  two-dimensional  spectrum  is 


$(2)(tc)  = 1.9xl05  k"4*95,  .0026  < K < .02  ft  1 


$(2) (k)  = . 53xl0~3  K“4,  .02  < k < 60  ft  1 


(5.4) 


yields  a 2 = 6.75  ft2  with  a peak  $(.02)  = 3300  ft4. 

2 

In  the  wavelength  region  1 to  50  cm,  recent  measure- 
ments^5 on  the  NELC  Tower  indicate  that  the  ocean  spectrum 

-a  15 

follows  a power  law  < . SRI  wire  gauge  measurements 

indicate  that  the  one-dimensional  spectrum  can  be  approxi- 
mated by  the  power  law 


${1)(k)  = k$(2)(k)  = 5xl0-3  k~ 3 1 -125  < K < 2u  cm-1 

(5.5) 

(1)  3 

where  $ (k)  is  the  one-dimensional  spectrum  (cm  /radian). 
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I I 
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230  22(f  2(0°  200° 


FINAL  SMOOTHED  SPECTRUM  FROM  SWOP 
( K,a)  FT.4 


K rad./ft. 

ONE-DIMENSIONAL  WAVENUMBER  SPECTRUM  TAKEN  ALONG  NORTH-SOUTH 
DIRECTION  IN  SWOP  CONTOURS  OF  FIGURE  5 

FIG.  6 
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RRI  measurements1  using  a TV  technique  indicate  an 
ocean  spectrum  approximated  by 


$(1>(k)  = 7. 82  xlO-5  <“4,  .125  < K < 2 tt  cm-1 


= 0,  otherwise 


(5.6) 


It  is  interesting  to  compare  the  rms  sea  height 

(o  ) and  rms  sea  slope  for  this  small  sea  roughness 

z ssa. 

regime  to  the  corresponding  SWOP  rough  sea  parameters. 

For  the  RRI  spectrum 


2ir 

az  = 2lT  f d 

.125 


d<(7. 82x10  5 <“4)  = 8.28x10  2 cm2 


a = *289  cm 
z 


(D  = .0621. 

wsea 


(5.7) 


A correlation  length  (S.)  appropriate  to  this  measurement 
is  estimated  assuming  a Gaussian  distribution  of  amplitudes, 
so  that 


®sea  * 


(5.8) 


yielding  a correlation  distance  of  9.21  cm.  For  the  SRI 
spectrum  o_  = 1 cm,  = .351  and  £ = 5.71  cm. 

Z 56a 
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In  Table  I are  collected  the  statistical  parameters 
for  the  four  ocean  spectra  discussed  and  for  three  radar 
wavelengths  of  interest,  UHF,  C band  and  X band.  This 
table  can  be  used  to  determine  which  of  the  radar  scatter- 
ing models  discussed  is  appropriate  for,  (i)  the  radar 
wavelength  being  used,  60  cm,  6 cm  or  1.6  cm,  and  {ii ) the 
sea  state  as  described  by  the  four  different  spectra.  It 
is  clear,  for  example,  that  only  the  scattering  model  in 
Appendix  E with  koz  >>  1 and  k£  >>  1 is  appropriate  for  the 
conditions  of  SWOP  for  any  of  the  radar  wavelengths . 

In  Figure  7 the  cross  section  per  unit  area  is  plotted 
vs.  angle  of  incidence  for  three  distinct  sea  states.  In 
the  SWOP  measurement,  as  pointed  out  above,  kc2  >•>  1 and 
k £ >>  1,  so  that  the  sea  state  is  heavy  and  the  scattering 
is  quasi-spectular.  From  Appendix  E,  the  cross  section  for 
arbitrary  radar  wavelength  is 


QS  - S 


4 

sec  Qc 
sea 


exp  - tan2©  / (D  2 

o sea 


(5.9) 


The  second  sea  state  considered  is  that  of  an  e^uil- 
briuro  or  Phillips  spectrum.  Perturbation  theory  for  hori- 
zontal and  vertical  polarization  from  Appendix  D yields 


— 


TABLE  I 


A=60  cm 
k=.1047  cm 


A=6  cm 
k=l .04?  _m" 


A*1 . 5 cm 


SWOP 

SRI 

PHILLIPS 

RRI 

o_  (cm) 

c* 

65.7 

1.002 

.448 

.289 

k0z 

6.88 

.1049 

.0469 

.0303 

® sea 

.0915 

.351 

.157 

.0621 

l (cm) 

1440 

5.71 

5.71 

9.31 

k£ 

150.8 

.598 

.548 

.975 

cz  (cm) 

65.7 

1.002 

.448 

.289 

koz 

68.8 

1.049 

.469 

.303 

®sea 

.0915 

.351 

.157 

.0621 

A (cm) 

1440 

5.71 

5.71 

9.31 

k£ 

1508 

5.98 

5.38 

9.75 

oz  (cm) 

65.7 

1.002 

.448 

.289 

k°Z 

275.3 

4.20 

1.88 

1.211 

® sea 

.0915 

.351 

.157 

.0621 

^ (cm) 

1440 

5.71 

5.71 

9.31 

kA 

6033 

23.9 

23.9 

39.0 

The  statistical  parameters,  a = rms  height  of  the  surface, 

* 5 correlation  distance , <Dsea  = rms  surface  slope,  and  k = 
radar  wavenumber  ( = 27r/A)  for  the  four  ocean  spectra,  SWOP,  SRI 
Phillips  and  RRI  for  three  radar  wavelengths  of  interest  UHF 
(60  cm),  C band  (6  cm),  and  X band  (1.5  cm)  are  listed. 
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S5T'  ~ Tl 


0°  = (pi)  = 16ir  k4  cos4©„  $(2k  sin©  , 0) 

HH  Vs  /HH  ° ° 


aW  E (s  J = 16lT  k4  U + sin20o]2  ${2k  sin0o'  0)  * 


(5.10) 

The  distinction  between  vertical  and  horizontal  polarization 
shows  up  clearly  for  large  angles  of  incidence  in  Figure  7 . 

The  final  example  is  the  physical  optics  approximation 
made  in  Section  III  for  a slightly  rough  sea.  Using  a Phillips 
spectrum  yields  the  cross  section 


= 16lr  k4  exP  I“(2ko  cos©  )2]  <M2k  sin©  , 0) 

ir  • U * o ZO  O 

(5.11) 

which  in  Figure  7 is  seen  to  correspond  closely  to  the  cross 


section  for  horizontally-polarized  radar  in  Eq.  (5.10). 
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VI.  CURRENT-INDUCED  PERTURBATION 

The  composite  model  discussed  in  Section  IV  enables 
one  to  separate  the  variable  surface  effects  from  the  direct 
scattering  effects  in  first;  order.  For  example,  consider 
the  surface  current  generated  by  an  internal  wave  propagating 
at  50  cm/sec.  A 64  cm, surface  wave  has  a group  velocity 

equal  to  the  internal  wave  phase  velocity.  A linear  analysis, 

6 7 

e.g. , Thomson  and  West  or  Case,  Watson  and  West 

would  predict  a resonant  transfer  of  energy  from  the  internal 
wave  to  the  surface  wave  spectrum.  The  waves  in  this  reson- 
ant region  would  therefore  be  strongly  perturbed.  The  visi- 
bility of  this  perturbation  to  short  wavelength  radar,  e.g., 
X-band  corresponding  to  3.6  cm  surface  waves,  would  be  given 
by  the  perturbation  in  the  long  wave  spectral  function  F($;0q) 
in  Eq.  (4.18) . 

In  addition  to  the  direct  variation  of  F(t;0Q)  (when 
the  long  waves  are  perturbed  by  a current)  is  a cascade  of 
the  perturbation  to  shorter  wavelengths.  In  this  process, 
the  current  modulates  the  longer  waves  which  in  turn  modu- 
late the  shorter  waves,  thereby  inducing  a variation  in 
4>(-Ak,r)  in  Eq.  (4.18). 

To  examine  the  effect  on  the  clutter  cross  section  of 

the  current-induced  spectral  perturbation,  we  follow  Watson 
18 

and  West  and  introduce  the  power  spectrum  of  the  surface 
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displacement  $(x,t)  at  a position  x 


F (x,  jc)  = | ^<a*(£+A/2)  a(k-A/2))>  exp  (iA«<) 


(6.1) 

where  a(k)  is  the  time  dependent  amplitude  of  the  complex 
surface  displacement  Z (x) , 


Z(x)  = 


. Z 


a(k)  exp(ik»x) 


(6.2) 


and 


C(x,t)  = | (z  (x)  - Z*  (x) ) . 


(6.3) 


To  use  the  continuum  wavenumber  variables,  we  define  the 
spectral  function,  over  a surface  area  S, 


->•  *>  Q -V  -> 

<i>(x,k)  = *•  F (x,k) 

(2tt  ; 


(6.4) 


with  the  normalization 


f d2k  $(x,k)  = <C2(x))> 


(6.5) 


If  the  mode  amplitude  for  the  complex  surface  dis- 
placement  in  the  presence  of  a surface  current  [a (k) ] is 
distinguished  from  that  in  the  absence  of  the  current  [A(k)], 
then  introducing  an  envelope  function  G£(x,t) , the  two  ampli 
tudes  are  related  by 


a (k)  = G £(x,t)  A (k)  . 


(6.6) 


If  the  surface  current  is  assumed  to  be  "turned  on"  at  time  t=0 
then  Gj£(x,o)  = 1.  We  can  employ  these  amplitudes  to  construct 
the  "ambient"  power  spectrum,  ie.,  the  power  spectrum  with  no 
current  present, 


Fa(x,1c)  = j £<A(M/2)  A (M/2)>  exp(i^*x) 
A 

The  perturbed  spectrum  is  discussed  in  Watson  and  West 


(6.7) 
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and  can  be  written  in  terms  of  the  continuous  spectrum  $(x,k) 
in  the  weak  current  case  in  the  following  forms: 

(i)  The  near  resonance  case; 


Ail ~ i(u  ict)  cos20  sinx?  + i(U  kt) 2 cosx£  (1 

i>  (h  2 o 2 ° 


tAp)Vsp 


+ 2 


.1/8sinxd 


2S  sin kE,  I /sinx£ 

Jr 


ft)2  nul  * 1 

UnV  [tKUo  sinKC] 


(6.8) 


(ii)  The  far  from  resonance  case; 


j,  U cosx£ 

5$ (x,k)  ~ o 


$ (k)  2 (CT-C  cos0) 

a JL  y 


jCT  (cos  0-2)  + 2 C cos©  (l-cos0)  + C cos20j 
II  g g J 


k 3^  ^ V 


UQ  COSX^  COS0 
, CT  - C COS0 

*■  i g 


(6-9) 
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and  $>  (L)  is  the  ambient  long  wavelength  spectrum. 

cl 

The  variation  in  the  clutter  cross  section  induced  by 
the  surface  current  can  be  obtained  from  Eq.  (4.18)  to  first 
order  as 
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resonant 


* non-resonant 


(6.13) 


The  first  term  in  Eq.  (6.13)  is  the  induced  perturbation  in  the 
long  wave  spectrum  by  the  surface  current.  Using  Appendix  B, 
we  may  approximate  this  term  by 


) 

1 - expf-  (K cU2(E'  - E )] 

' resonant 

w Cm  Cm  Cm  J 

(6.14) 


so  that  using  Eqs.  (4.13)  and  (6.7)  we  obtain 


(£) 

resonant 


(6.15) 


Equation  (6.15)  may  be  evaluated  by  substituting  for  the  resonant 
perturbation  of  the  long  wave  spectrum,  Eq.  (6.8). 

As  an  application,  we  choose  for  the  ambient  spectrum, 

. . 13 

Phillips  equilibrium  spectrum. 


V*>  = I 


k within  90°  of  wind  direction 


B = 4x10 


(6.16) 


Thus , 


* =r- fc"  * <wl 


-4  cos© 


(6.17) 


II 


X8 

Following  Watson  and  West  we  choose  parameters 
to  correspond  to  typical  oceanographic  conditions: 


_2 

UQ  = 10  m/sec 


k = 10  2 m-1 


k = 8.66  m 
r 


-1 


Cj  = .50  m/sec 


0 = 20' 


(6.18) 


I 


Here  kr  is  resonant  wavenumber  defined  by  C1  = cos0.  The 
strength  parameter  Sp  is 


i a 


S = 4.609x10 
P 


(6.19) 


the  resonant  time  is 


? 

x = r — = 23.09  sec 
r KrUo 


(6.20) 


and  the  width  of  the  resonance  region  is. 


1/2  -1 
Ak  = kS  = 2.147  m 
P 


(6.21) 


40 
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For  the  Phillips'  spectrum  the  term  involving  an  integral 


over  L in  Eq.  (6.10)  is  negligible  compared  with  the  other  two. 

3(£n«  ) 

In  Eqs.  (6.8)  - (6.10)  the  term  involving  k — appears  to 

dominate  ^except  for  t >>  Tr  in  Eq.  ( 6 . 8 ) J . This  was  pointed 
18 

out  by  Milder  in  his  discussion  of  spectral  perturbations. 

For  a spectrum  of  the  form  the  relative  contribution  of  the 
first  and  last  terms  in  Eq.  (6.10)  is  approximately  ^ , which 
in  two  dimensions  is  .125. 

For  the  parameters  given  by  Eqs.  (6.18)  - (6.21),  Eq. 
(6.15),  for  an  adverse  current  becomes. 


resonant  = 1 ' ) 

= 1 - exp^-4k^az2cos20o^l-.Ol^p-^simc^  ] (6.22) 

The  nonresonant  contribution  to  the  variation  in  the  clutter 
cross  section  is  given  by  Eq.  (6.10)  with  k >>  kr  ~ L 


nonresonant 


-0.06  coskx 


(6.23) 


The  (&n<I>a)  term  provides  the  dominant  contribution  to  the 

x 

modulation.  The  fractional  change  in  the  clutter  cross  section 
using  Eqs.  (6.22)  and  (6.23)  is 


.06 


coskx  - exp 


► 


4 (ka  )2cos2 
z 


e4- 


(.01) 


We  see  from  Table  I that  ko  for  X-band  radar  is  1.88 

z 


for  a Phillips  spectrum,  so  that  the  resonant  contribution  is 
maximum  for  0Q  - 0°  and/or  early  times.  The  fractional  change 
in  cross  section  for  small  angles  of  incidence  is 


i 

f 


6o 

a 


1 


.06  COSKX 


(6.25) 


so  that  the  total  cross  section  varies  by  a factor  of  two  at 
X-band  frequencies. 

C-band  radar,  however,  has  ko  = .0469  for  a Phillips 

z 

spectrum  so  that  for  the  same  perturbing  current 

- -.06  coskx  (6.26) 


and  the  variation  in  clutter  cross  section  is  12%  peak  to  peak. 

Another  source  of  variation  in  the  clutter  cross  section 
is  that  due  to  changes  in  angle  of  incidence  of  the  radar 
(ie.,  changes  in  surface  slope).  Labeling  the  fractional 
change  in  clutter  cross  section  in  Eq.  (6.24)  as  being  current 
induced,  we  write  the  total  variation  as 


a V o / 

' current 

induced 


+ 


rms 

slope 


(6.27) 


From  Eqs.  (4.15),  (4.16)  and  Appendix  B we  have  that  the  depend- 
ence of  the  long  wave  scattering  on  the  angle  of  incidence  is 
to  lowest  order  in  the  rms  slope 
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Assuming  a modulation  function  M(k,x)  such  that 
6$ (k,x)  = $ (k)  M (k ,x) 

— •—  a 


(6  32) 


eq.  (6.30)  becomes 


rms 

slope 


-1.7x10 


sinic£ 

e2 

xnc 


(6.34) 


or  in  terms  of  the  current  itself 


(?)  = • 
rms  r inc 

slope 


The  relative  perturbation  in  the  clutter  cross  section,  therefore, 
2 

falls  off  as  0.  and  is  proportional  to  the  perturbing  current. 
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VII.  DISCUSSION  AND  CONCLUSIONS 

The  ocean  surface  has  at  least  two  classes  of  roughness: 

one  class  where  ko  >1  and  one  where  kc  <1.  In  the  NELC 

2 z 

Tower  measurements,  six  radar  frequencies  (P,  L,  S,  C,  X,  K) 
were  used.  For  the  sample  scattering  calculations,  the  nominal 
values  X = 60  cm  (P) , X = 6 cm  (C)  and  X = 1.5  cm  (X)  were 
used.  We  assume  that  the  ocean  spectrum  is  a composite  of  the 
rough  ocean  (SWOP)  and  the  slightly  rough  ocean  encountered  in 
the  NELC  measurements.  If  these  types  of  roughness  are  present, 
they  result  in  scattering  having  characteristics  of  each.  To 
a good  approximation,  one  can  add  the  incoherent  cross-sections 
directly  to  obtain  the  total  incoherent  cross-section  if  the 
processes  are  statistically  independent.  The  large  scale 
roughness  produces  strong  components  near  normal  incidence 
which  are  identical  for  vertical  and  horizontal  polarization. 
This  is  the  so-called  quasi-specular  component.  Small  scale 
roughness  contributes  to  the  tail  (0Q  near  90°)  and  is  referred 
to  as  the  diffuse  component. 

In  Table  I the  statistical  parameters  for  four  ocean 
spectra  and  three  radar  wavelengths,  including  UHF,  C banc  and 
X band  are  listed.  From  this  Table  it  is  clear  that  the  appli- 
cability of  a specific  model,  ie.,  set  of  assumptions  for  the 
radar  scattering,  is  dependent  on  both  the  sea  state  and  radar 
wavelength.  In  Table  II  the  scattering  models  considered  in 
this  paper  are  listed  where  most  appropriate  for  a particular 
sea  state  and  radar  wavelength. 
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TABLE  II 


From  Section  VI  it  is  clear  that  the  "visibility"  of  a 
surface  wave  spectral  perturbation  is  dependent  on  the  sea 
state  and  therefore  on  the  radar  frequency  with  which  one 
views  the  sea  surface.  In  the  example  considered  the  visibility 
of  the  current  induced  perturbation  was  reduced  by  an  order  of 
magnitude  when  one  changed  from  X to  C band  radar  in  a Phillip's 
saturated  sea. 

The  variation  in  the  clutter  cross  section  in  lowest 
order  consists  of  three  terms. 


/SF\ 
\F  ) 


resonant 


nonresonant 


(¥) 


rms. 

slope 


The  first  term  is  the  direct  mechanical  interaction  of  sea 
waves  resonant  with  the  surface  current.  The  second  term  is 
the  perturbation  in  the  short  surface  waves  due  to  both  the 
nonresonant  interaction  with  the  current  and  wave-wave  inter- 
actions of  the  scatterers  with  the  sea  waves.  The  final  term 
is  due  to  the  variation  in  the  angle  between  the  sea  waves  and 
the  incident  radar  and  decreases  as  the  inverse  of  the  square 
of  the  angle  of  incidence. 


Biilb&to 


V-' 


P4, 
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(A3) 


This  defines  a parabola  within  the  geometric  shadow  which 
is  illuminated  due  to  diffraction. 

The  maximum  area  illuminated  due  to  diffraction  would 
be  twice  the  wave  height,  located  at  the  wave  trough,  i.e.. 
Ax  = X /2,  therefore, 

SGcl 
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r 


2H  = = /ArXsea^ 


(A4) 


This  equation  establishes  criteria  for  the  threshold  of 
shadowing,  that  is, 


8H  VXr  > Xsea; 


8H  /Xr  >>  Xsea; 


illuminated 
(no  shadow) 


geometric  optics 
(no  diffraction) 


(A5) 


In  terms  of  oceanographic  parameters,  such  as  the  slope 

0 = 2-irH/A  , we  can  write 

sea  ' sea 


VXsea  < ® sea  * k ® 


2 

sea 


(A6) 


where  the  rms  slope  0 has  replaced  that  for  a single 


sea 

wave  0 . For  a saturated  sea 

sea 


® L.  = -°06  iog 


sea 


/ jmax\ 
\kmin/ 


(A7) 


so  that  reasonable  numbers  for  the  mean  square  slope  are 
like  .025  to  .05.  Substituting  these  values  into  Eq.  (A6) 
yields  .005  £ A /A  .01,  so  that  for  X-band  radar 

(Ar  ~ 3 cm)  the  sea  waves  which  give  shadowing  are  ^ 3 to 
meters.  These  criteria  establish  the  region  of  the  wave- 
number  spectrum,  below  which  shadowing  can  be  neglected 


I*,  jnwniMKi 


Appendix  B - Averages  of  Statistically-Separable  Surfaces 


Consider  the  average  of  the  series  expansion  of  the 
characteristic  function, 

2 

^exp  [iK(Z-Z')]^  = 1 + iK<^Z-Z'^  - |-<^(Z-Z')2^ 

- 1 r <(z-z')6>  + Si  <tz-z')4>+  . 


The  odd  moments  in  B1  vanish  for  a Gaussian  distribution 


<<Z-Z')2>=  2az2  [1  - Ez(p-p1)] 


<<Z-Z')4>=  12az4  [1  - £ z (p-p* ) ] 2 


so  that  the  series  expansion  becomes 


<exp[iK(Z-z')]}>  = 1 - (Kaz)2  (1-Ez)  + 


M 


r-  «-*.>*  + 


It  is  clear  that  the  series  in  B2  can  be  summed  to  yield 
the  exact  result 


^exp  [iK(Z-z')]^=  exp  [-(Kaz)2(l  -E^jJ 


51 


Similar  exact  calculations  may  be  made  to  obtain  the 
following  results; 


<(zx  + Z'x, ) exp[iK(Z-Z')J^  = 


= -i2K0 


2 Ih 

z 3x 


exp^-(Kaz)2(l-Ez)] 


(B4) 


4Zx  + Zx,)  expjiK(Z-Z#)J^  = 2 
<ZXZ^,  (Zx  + Z',)  exp  j^iK  ( Z -Z  1 )] ^ 


4 -^z)24^ 

X 


(?3r)  ^p[-(K02)2(l-£z)j 


(B5) 


* i2Ka 


2 

z 3x 


<Z2  Z^?  exp[iK(Z-z')]> 


exp  [-(Kaz)  2 (1“  E )] 
(B6) 


, 4 32E  v 2 

4_  43  z I 


+ 4(ko  ) a _ 
z 2 3x2 


-2  (Ko  )2a2 
z Z, 


fr) 


<^zvz^f  exp[iK(Z-Z')J^>  = - 


x exp  (K  az  ) 2 ( 1-  Ez  )J 


°z2  -f-  + (Koz)2°z2 
2 dxZ  z 2 


(B7) 


exp[-(Koz)2(l-E2)] 


(B8) 


I 

i 
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C.  PHYSICAL  OPTICS  SCATTERING  FOR  A ROUGH  SEA 

In  Ref.  (8)  the  physical  optics  approximation  was 
used  to  evaluate  the  Stratton-Chu  scattering  integral. 

This  tangent-plane  approximation  requires  the  radius  of 
curvature  at  all  points  of  the  sea  surface  to  be  greater  than 
the  radar  wavelength  Ar.  No  restriction  was  made  concerning 
the  mean  square  roughness  height  compared  to  Ar>  in  contract 
to  Ref.  (1)  where  H/xr  <<  1.  Otherwise  the  model  is 
similar  to  that  described  in  the  previous  sections,  i.e., 
any  correlation  lengths  are  much  greater  than  the  wave- 
lengths of  the  scattering  centers.  For  a Gaussian  correla- 
tion function,  it  is  possible  to  evaluate  the  incoherent 
scattering  cross  section  exactly.  The  resultant  backscatter 
cross  section  per  unit  area  is,  using  l as  a correlation 
length, 

•jjr  = (k£)  2 sec20Q  4»  (ct,  B) 
a = (2kcr_  cos0  ) 2 

£*  U 

8 = (k l sin0Q) 2 

00 

-7-  exp(— B/n) . (Cl) 

n=l  n*n 

Figure  8 depicts  the  incoherent  scattering  function  ^(a,g) 
extracted  from  Ref.  (8). 
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FIG.  8 INCOHERENT  SCATTERING  FUNCTION  FOR  A SURFACE  WITH  A GAUSSIAN  CORRELATION 
FUNCTION 


D.  PERTURBATION  THEORY  FOR  SCATTERING  FROM  A SLIGHTLY 
ROUGH  OCEAN 

In  previous- sections,  the  physical  optics  approxi- 
mation was  made  which  is  valid  when  the  radius  of  curvature 
of  the  ocean  surface  is  large  compared  to  the  radar  wave 
length.  In  this  section  we  relax  this  assumption  and 
require  that  the  roughness  height  be  small  (k a_  <<  1) 
and  that  the  surface  slopes  be  relatively  small  (V£  < 1) . 

A smoothly-curving  surface  with  small  slope  will  satisfy 
the  physical  optics  radius  of  curvature  requirement  for  the 
high-frequency  limit  but  eventually  this  condition  will 
fail  a.  the  frequency  is  decreased. 

Consider  the  scattering  geometry  shown  in  Figure  9 . 
For  small  surface  fluctuations  (S)  about  the  reference  ocean 
surface  {SQ) , we  expect  small  perturbations  in  the  scattered 
field.  When  Sq  is  illuminated  by  a plane  wave 


. ^ A A A ^ 

E^(r)  = (xA1+yA2+zA3)  exp(-i)c-r)  , Jc  = (k1,k2,k3) 


then  the  scattered  field  is  (with  SQ  large  and  infinite 
ocean  conductivity) 


£g(r)  = (-xA1-yA2+zA3)  exp  J-i  (k^x+k^-k^)  . (D2) 
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SCATTERING  GEOMETRY  WITH  INCIDENT  WAVE  IN  THE 
(X,Z)  PLANE  AND  A SMOOTH  OCEAN  SURFACE  (SQ) 
AND  DISPLACED  SURFACE  (S) 


This  solution  satisfies  the  tangential  boundary  conditions 


z * £ = 0 on  SQ,  £ = i?  + i?s 


■z  x $ (?•)  = z x £ (?•>  - *tj;o) 

so  o o m 


?'  c 
o o 


We  can  also  view  the  above  scattering  problem  in 

the  following  way.  The  incident  field  induces  a magnetic 

current  iti°>  on  the  surface  S which  qives  rise  to  a field 
m o 3 

above  the  surface.  The  total  scattered  field  is  obtained 
by  summing  all  such  magnetic  surface  currents  over  SQ. 

This  is  facilitated  by  using  the  magnetic  dyadic  Green's 
function  for  a half  plane 


-2V  x 


When  the  surface  is  rough  the  scattered  field  will  be  per- 
turbed from  that  given  by  Eq.  (D4)  . We  calculate  the  per- 
turbed field  by  finding  the  equivalent  surface  currents  on 
S which  account  for  the  current  on  S at  the  height 


C = C(x',y')  = 


The  currents  on  SQ  are  found  by  using  a perturbation 
expansion  of  the  scattered  field  in  powers  of  the  surface 


Wil&riturtdMd! Hill .mini Hill. I ... 


height  C-  This  expansion  satisfies  the  boundary  condition 


n(r* ) x [^s (r ' ) + = 0 


Let  the  scattered  field  be  represented  by  the  series 


o c 

gs(?,)  = X! g{n)  (?,)  ; |fi( 


(n)  . . ,n 

~ u 


where  each  e(n)  satisfies  Maxwell's  equations.  In  addition, 
each  term  will  be  expanded  in  a Taylor  series  about  r^. 

To  first  order  the  scattered  field  is 


Es(?')  » £(o)  (?^)  + c ~t  2(0)  (P)  + fi(1)  (?;) 

£(o)  (r')  = X exp(ilt.r')  » Xd-ik,*;)  exp (-iJt.?^)  . (D6) 


With  the  equation  for  the  free  surface  z*  — ^(x',y')  = 0, 
the  unit  vector  normal  to  this  surface  is  to  first  order 

id’)  - a a.. 3..  (D7, 

|v(2*-0|  I z~xcz ?~y^y » j 3X  y 

When  Eqs.  (D6)  and  (D7)are  used  in  the  tangential  boundary 
conditions 


;[i  exp(-k-r^)  +E(o)  (r^)j  = 0 
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W'l'iqiHiiHil 


we  obtain 


zx  -ik3cA  exp(-it-r^)  + ? g|r  E(o) (r^)  + E 


(1) 


(iCjj.+JCy.)  * exp(-iic-r^)  + l(0)  (r^) 


(D8) 


From  Eq.  (D8)  we  readily  calculate 


z x 


£{1)  (r^)  = 2^i 


A A 


A A 


:•'] 


|ik3i;(yAl_xA2)  + A3{x^y,_y?x,)  I exP<-ik*ro} 

(D9) 

The  first  order  magnetic  current  is  then 


^m  )(r^  2 ik3^A2+A3^y ' ' + yf^k3^Al+A3^x*  ^ j exp(-k-r 


')  (DIO) 
o 


which  could  be  used  in  Eq.  (D4)  to  calculate  2^  (f) . 

However,  it  will  be  more  convenient  to  calculate  the  magnetic 
field  using 


5s{?)  = SiT  Jds^T+  pj  G(r-r^)  . (xx*yy)  -^(r^ 


(Dll) 


Using  the  far  field  approximation  to  simplify  the 
Green's  function  in  Eq.  (Dll)  along  with  Eq.  (DlO),  we  obtain 
for  a ~ectangular  surface  (SQ  = xQyo) 


a(o)(?>  = exEj-ii.r) 

2tt1(j3vi  r 


where 


inc(qi  -#)(T- 


A a\  A A 

rrj-  ( iA2~ j A3) 


(D12) 
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q = xqx  + yq2  , = k-j-kr.x  , q2  = k2-kr.y  . (D 


Using  the  spectral  representation  for  the  random  surface 
displacement 


S(r’)  = 


d$(<)  exp (-iK’ r^)  , £ = X|C1+^K2 


and  the  vector  §(<)  defined  by 


§(ic)  = acB1  + yB2  = x(k3A2+A3<2)  - y (k3A1+A3ic1)  (D14) 


the  next  order  term  in  the  expansion  of  the  magnetic  field  is 


3(1)(J)  = ^ exp^-ikr)  jcs  exp  (iki-r ' j/d*  (k)S  '?.) 

X exp (-1 (k+k) *r^)  So 

ar  ing  cut  the  surface  integration  in  (D15) , we  find 
for  the  integral  terms,  call  it  the  vecuor  5(r), 


3(r)  = xQy 


Oo 

r ^ ^ x y 1 

0 / d$(ic)  §(k)  sine  (q^-i^)  ~y  sine  (q2~i<2)  -Tpj 


(D16) 


3{i)(?)  = exH^ikr),  * 


(D17) 


The  radar  cross  section  in  terms  of  the  magnetic 


field  is,  to  first  order. 


o = 4 -nr 


2 4,r2  f ,M°)  ,2  . ✓ ,*•<1)  ,2S 

1 <'  1 > 


where  - 0 and 


H 2 , £ E2  . 


(D17) 


O y 

The  first  order  contribution  to  the  cross  section  is 


r{1)  = 

ttE  * 


(<|ff|2>  - <i?-«i2>) 


(D19) 


where 


00 

<J0|k>  = (x0yQ)2  J <d$(£)  d$*(£-)>  $(£>.&(£') 


x sine 


3 [te]."*!)  x]  sinc  [(q2"lc2)  x]  ' (D20) 


For  a spatially-homogeneous  spectrum  of  ocean  surface  waves 

<(d$(<)  d$*  <£*)/>  = d<  dje'  $(<)  5(k-k*)  . (021) 

When  Eq.  (D21)  is  used  in  Eq.  (D20)  and  the  ocean  spectrum 
does  not  change  too  rapidly  about  q,  then  the  approximate 
value  of  integral  (D20)  is 
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<|D|2^  - |f(q)  |2  $(q)  (2  ir)  2 xQye 


(D22) 


Similarly, 


<|r.S|2>  m J2  «($)  (2ir) 2 xoyQ  . 


(D23) 


o(1)  then  becomes 


(1)  4k  ,,  v 
cr  * — x-  (2  tr) 


2 x0y0  |fl(|)  |2-  |r.S(4)  |2  ]«(q)_ 


(D24J 


We  now  specialize  to  backscatter  (r  » -k)  and  hori- 
zontal polarization  with 


A1  ’ A3  ’ 0 


* -k  sineo,  kj  «»  0,  k^  * -k  cosSj. 


From  Eq.  (D13) 


q^  * k^  - kr*x  * -2k  sine< 

A A 

q2  * k2  - kr.y  ■ 0 


and  Eq.  (D14) 
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B-^q)  = k3A2  + A2q2  = -k  cos 0qEq 


B2(q)  = k3Al  + A3ql  = 


= 0 . 


Substituting  the  above  conaitions  for  backscattering  in 
Eq.  (56)  yields  the  first  order  incoherent  cross  section 
per  unit  area  for  horizontal  polarization 


4 4 

^ - — = 16 tr  k cos  eQ  $(-2k  sin0o,O) 
o^o 


(D25) 


Next  we  specialize  to  backscatter  and  vertical  polarization 

, A A 

for  which  EQ  = xA^^  + zA^  and  A2  = 0,  all  other  quantities 
being  the  same  as  in  the  preceding  example.  With  these 
values  the  incoherent  cross  section  per  unit  area  for 
vertical  polarization  is 


° — = 16tt  k4  (1  + sin2eo)2$(-2k  sineQ,0) 
o^o  ' "* 


(D26) 


Using  Eq,  (D12)  for  S ° , we  determine  the  zero**  order, 
or  coherent  cross  section  for  both  vertical  and  horizontal 
polarization,  to  be 


o'  y 2 2 

— - = (kxoyQ)  cos  0q  sine  (kxQ  sin0Q) 
o ■ o 


(D27) 
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E.  QUASI -SPECULAR  SCATTERING  FROM  VERY  ROUGH  SEA 

When  surface  roughness  is  large  compared  to  radar 
wavelength  (ka_  >>  1) , the  scattered  power  becomes  mor<; 
incoherent.  For  kr?  > 5,  kr  being  the  radar  wavenumber  2ir/\r, 
the  coherent  power  is  negligible.  The  following  analysis 
assumes  that  all  scattering  comes  from  specular  reflections 
so  that  the  backscattered  power  is  perfectly  reflected  and 
consequently  there  is  no  cross  polarization. 

When  surface  roughness  is  isotropic,  the  correlation 
length  for  surface  waves  is  much  less  than  either  dimension 
of  the  illuminated  area.  Neglecting  multiple  scattering 
and  shadowing  allows  for  stationary-phase  integration  of 
the  surface  scattering  integral.  This  integration  yields 
the  incoherent  backscattering  cross  section  per  unit  area: 


4 

sec  9 


-^-2  exp(-tan2eo/®2ea) 

eoa 


where 


2 f2°z\ 

® sea  2 VT/  is  the  raean  square  total  slope  of 
the  rough  surface 

0 -r2/!2  r: 

* 2 correlation  length  for  p( r)  = e 7 S 1 - ~ 


r = [(xfx2)  + (y1~y2)  I*5 


for  a Gaussian  surface-height  joint  probability  density. 
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